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Abstract 

We calculate the maximum mass of the class of compact stars described by Vaidya- 
Tikekar [24] model. The model permits a simple method of systematically fixing bounds 
on the maximum possible mass of cold compact stars with a given value of radius or 
central density or surface density. The relevant equations of state are also determined. 
Although simple, the model is capable of describing the general features of the recently 
observed very compact stars. For the calculation, no prior knowledge of the equation 
of state (EOS) is required. This is in contrast to the earlier calculations for maximum 
mass which were done by choosing first the relevant EOSs and using those to solve the 
TOV equation with appropriate boundary conditions. The bounds obtained by us are 
comparable and, in some cases, more restrictive than the earlier results. 
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1 Introduction 

In relativistic astrophysics the problem of maximum allowed mass of a very compact star 
is an important issue. The upper bound on mass, in the case of an ultra-compact object, is 
crucial for making the distinction between a black hole and a 'normal' star. However, the 
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estimation of limiting values of masses is very much model dependent. The maximum mass- 
radius calculation is essentially based on the assumed equation of state(EOS), although 
some physical processes which can lead to the instability of the star should also be taken into 
account. Majority of calculations, based on different neutron star EOS, put a constraint 
on maximum allowed mass within the range M ~ (1.46 — 2.48) M & and radius b ~ (9 — 
11.7) km (see Table 2 of ref. [15]). However, one can show that above a certain density, 
it is possible to put a constraint on the maximum allowed mass irrespective of the exotic 
nature of the EOS used to describe such stars [1]. For a non-rotating neutron star, Rhoades 
and Ruffini [2] have shown that this maximum mass can not be larger than ~ 3.2 M . 
However, data obtained from new generation satellites like 'Chandra' predict much lower 
masses for pulsars than predicted theoretically in the framework of neutron star EOS. 
Thorsett and Chakrabarty [3] have pointed out that masses of the majority of the pulsars 
lie within a narrow range of (1.35 ± 0.04) M which also needs an explanation. The new 
class of compact stars may confirm the conjecture made independently by Witten [4] and 
Farhi and Jaffe [5] that quark matter (and not 56 Fe) is the true ground state of hadrons. 
Recent investigations show that a number of pulsars like Her X-l [6], 4U 1820-30 [7], [8], 
SAX J 1808.4-3658 [9], 4U 1728-34 [10], PSR 0943+10 [11], [12] and RX J185635-3754 [9], 
[13], earlier thought to be neutron stars, are actually good strange star candidates. Kapoor 
and Shukre [14] claim that none of the modern neutron star EOS is good enough to explain 
the observed mass-radius constraints of known pulsars. 

For more compact (strange) stars, one may make use of the MIT Bag model. The EOS 
for strange stars made of it, d and s quarks, can be written as [4], [5] 

P =l(p-4B), (1) 

where p is the energy-density, p the pressure, B the Bag constant. In the framework of the 

phenomenological bag model, Banerjee et al [19] have obtained a kind of Chandrasekhar 

limit on maximum mass for strange stars from a consideration of the binding energy. The 

maximum mass obtained in their calculation shows a scaling with respect to the Bag 

constant B. Recently, Dey et al [8] have obtained new sets of EOSs for strange matter 

based on a model of interquark potential which has the following features: (a) asymptotic 

freedom, (b) confinement at zero baryon density and deconfinement at high baryon density, 
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(c) chiral symmetry restoration and (d) gives stable uncharged /3-stable strange matter. 
These EOSs have later been approximated to a linear form by Gondek-Rosihska et al [20] 
which has the from, 

p = a(p-p ), (2) 

where po denotes the energy density at zero pressure and a is a constant. Making use of 
equation (2) Harko and Cheng [21] and Zdunik [23] have shown that for po = 4B (B = 
56 MeV fm~ 3 ), the maximum mass of a strange star is M max = 1.83 M . Interestingly 
the maximum mass and radius of a strange star is very much similar to that of a neutron 
star, the difference being in their mass-radius curves [16], [17], [18]. 

Naturally one may ask if it is possible to put constraints on the maximum allowed masses 
of ultra-compact stars without going into the details of the complex EOSs governing the 
configurations of such stars. The purpose of this paper is to investigate this problem 
for a class of stars based on a model given by Vaidya and Tikekar [24]. Exact solutions 
of Einstein's field equations, compiled by Finch and Skea [25] and Delgaty and Lake [26], 
reveal that although we have a large abundance of exact solutions (more than a thousand!), 
most of the solutions are not suitable for the description of realistic stars. In some earlier 
works [27], [28], it has been shown that the general solution obtained by Mukherjee et al 
[29] for the class of stars described by Vaidya-Tikekar [24] model, is very much relevant for 
the description of superdense stars. Vaidya-Tikekar [24] model has already been studied 
by many workers (e.g., [30], [31], [32], [33], [34]) for describing ultra-compact stars. In 
particular, Knutsen [34] for a specific value of a parameter (e.g., A = 2 in equation (5)), 
has shown that the maximum mass of a relativistic fluid sphere can not exceed 3 Mq 
for Baym et al's [35] EOS. Using the solution obtained by Mukherjee et al [29], which is 
valid for all physical values of the parameter A, it is possible to calculate analytically the 
maximum possible masses for these class of stars, given the necessary inputs. The model, 
though given in terms of simple analytic expressions, has been found to be able to describe 
the essential features of realistic compact stars. In particular, we have shown that the 
model can be used to describe pulsars like Her X-l and SAX J 1808.4-3658 [27], [28]. 

The presentation of the paper is as follows. In section 2, we give an outline of the model 
which will be relevant for the analysis of our work. In section 3, we obtain bounds on 
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the maximum mass for given values of radius or central density or surface density of the 
compact star. The results obtained are discussed and compared with earlier results in the 
last section. 



2 Vaidya-Tikekar model 



We consider a static, spherically symmetric star whose interior metric is given by 

ds 2 = _ e 2y(r) dt 2 + e 2»{r) dr 2 + ^2 + ^2 Qd(j) 2y (3) 

The matter within the star is assumed to be of perfect fluid nature and consequently we 
choose the energy-momentum tensor in the form 



Tij = (p + p)uiUj +pgij, 

where u l is the 4-velocity of the fluid, p is the energy-density and p the pressure. 
Starting with the ansatz, given by Vaidya and Tikekar [24], 



(4) 



a 2 M 



1 + Xr 2 /R 2 
1 - r 2 /R 2 



(5) 



Mukherjee et al [29] made use of the pressure isotropy condition to obtain the other 
metric function, 



=7 _ 



A 



cos[(n + 1)C + 5] cos[(n - 1)C + 5] 



n + l 



n — 1 



(6) 



In equation (6) Q = cos 1 z, where, z 2 = ^xx(l — jip) and n 2 = A + 2. 
The energy-density p and pressure p, in this model, are given by 



1 



R 2 (l - z 2 ) 
1 



1 + 



p = — 



(A + l)(l-z2) 

2zip z 



1 + 



(A + 1)V 



(7) 



(8) 



R 2 {1 - z 2 ) 

where ip z denotes differentiation of if) with respect to z. The total mass of the star for a 
radius b is given by 

(1 + A)6 3 /,R 2 



M(b) 



2(1 + \b 2 /R 2 Y 
4 



(9) 



The model has four parameters, A, 5, R and A. The conditions p = and 7 = —p, 
at the boundary, r = b, will determine two of these parameters (say A and 5). We need 
physical inputs to determine the rest. Assuming that A specifies the EOS for the given star 
in this model, any one of the following, e.g., the radius (6), central density (po) or surface 
density (pi,) will determine the parameter R. Thus the model gives a complete description 
of the star, including the radial dependence of density, eqn.(7) and pressure, eqn.(8), and 
determines the EOS. 

3 Relativist ic bound on mass 

A relativistic star, in general, satisfy the following conditions [1]: 

(i) Inside the fluid sphere the velocity of sound < v s (= \J~^) < 1- Since we are consid- 
ering only cold stars, the identification of the velocity of sound with (i. e . ; i n principle 
at zero temperature) is unambiguous. 

(ii) The fluid sphere should be dynamically stable against small radial perturbations. 



In addition, for a realistic star, the interior solutions must satisfy certain regularity 
conditions, namely, 

(iii) Both energy-density and pressure should be positive in the interior of the star. 

(iv) Pressure should vanish at some finite distance from the centre of the star, determining 
the radius of the star. 

We will make use of these principles to find bounds on the stellar configurations in 
Vaidya-Tikekar model [24]. From equations (7) & (8), we obtain 

dp z(i- rTO.A/O 2 - (i- z 2 ){i> z /i>) 



dp 



z(l-z 2 )(l + A)+4z 



(10) 



Condition (i), then, implies 

1 



l-z 1 



2z 



< — < 



tl> ~ l-z 2 



— +D 

2z 



(11) 



where, 



D 



4 + ^ + (l + A)(l-z 2 ) 



1/2 



Imposition of the condition (hi) demand that A > — 1 and 



*<-<i±*>. (12) 



Imposition of the condition (iv) gives 

Mzb) _ (l + A) 
ip(zb) 2z b 

where, 



(13) 



In ref.[29], combining the two constraints given by equations (11) & (12) an effective 
bound on tp z /tp was obtained as 



z 2 



1-D 

2z 



< t± < S1±A. (14) 

- V ~ 2z V ; 



Equating the two sides 

(a) at r = gives A > 3/17 and at 

(b) r = b yields 

b 2 A 2 + 5A + 12 - (17A 2 + 82A + 129) 1 / 2 

i^- 1 M5 + ^) ' ( } 

A simple calculation will, however, show that equation (15) is a weak bound and it 

does not ensure the restriction dp/ dp < 1 everywhere. For example, if A = 2, the above 

condition implies that 6 2 /i? 2 < 0.5. However, it can be shown that if 6 2 /ii 2 > 0.4235, 

the causality condition is violated. Therefore, a more stringent condition, necessary and 

sufficient to satisfy the causality condition should be given. It may be pointed out that 

out of 127 solutions of Einstein's field equations compiled by Delgaty and Lake [26] only 

16 qualify all the tests discussed in the earlier section. Out of 16, only 9 solutions have 

sound speed which monotonically decreases with radius, a criterion we intend to adopt. 

• We assume that the value of ^ is maximum at the centre. The condition ^ < 1 at 
the centre, then, yields 
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(j) Z0 > v / ATl±v / 2lATl , (16) 



where, 

°"(A + 1)- 

Now for A > 3/17, positivity of pressure (see equation (8)) demands that ^ must 
be a negative quantity which is possible if we choose only the negative sign in the 
expression on the right hand side of equation (16), i.e., 



i>z , ^ (i + A) 



V' 2 °- 2^A 
Again, from equation (6), we have 



v/XTT - V21A + 1 



(17) 



Y> 2 (n 2 - 1) 



sin[(n - 1)C + S] - sin[(n + 1)C + 5] 



(n + 1) cos[(n - 1)C + 5] - (n - 1) cos[(n + 1)£ + 5] 



(18) 



Combining equations (17) (with equal sign) and (18) at the centre, one can, thus, 
determine the limiting value of 6 for a given value of A. 

• Corresponding to the limiting value of 5, equation (13) can be used to calculate the 
maximal value of b 2 /R 2 for the same value of A. 

• From equation (9) the compactness of a star in this model is given by 

M(fc) = ^+A) 
b 2(A + i)' 

where, y = b 2 /R 2 . Clearly, the maximum value of y corresponds to the maximum 
compactness for a given value of A. 

The method employed here to calculate the maximum mass can be summarised as 
follows: 

(1) We specify a value of A (note that the same A may give different EOS in different 
stars). 

(2) We assume that the maximum mass corresponds to (^)o = 1> i-e., ^ = 1 at r = 0. 



This determines y m ax = (^)max and hence the maximum compactness u max = (x) 



max — V b ) max ' 



A 


( & \ 

Jmax 


( M\ 

\ 5 frnax 


M max /M & 








b = 10 km 


6 = 8 km 


6 = 6 km 


1 


0.4618 


0.3159 


2.14 


1.71 


1.28 


2 


0.4234 


0.3438 


2.33 


1.86 


1.39 


3 


0.3727 


0.3519 


2.38 


1.90 


1.43 


4 


0.3297 


0.3554 


2.41 


1.92 


1.44 


5 


0.2944 


0.3573 


2.42 


1.93 


1.45 


7 


0.2417 


0.3591 


2.43 


1.94 


1.46 


10 


0.1898 


0.3602 


2.44 


1.95 


1.46 


20 


0.1102 


0.3611 


2.44 


1.95 


1.46 


50 


0.0486 


0.3614 


2.45 


1.96 


1.47 


100 


0.0252 


0.3615 


2.45 


1.96 


1.47 


200 


0.0128 


0.3615 


2.45 


1.96 


1.47 



Table 1: Maximum mass (M max ) of a star for different radii 6 = 10,8, and 6 km and for 
different choices of the parameter A. 



(3) There is still one free parameter R. It can be specified by giving the values of one 
of the three: (i) radius 6, (ii) central density po or (hi) surface density pb- The relevant 
relations, in this model, are: 

3(A + 1) (i + A )(3 + Aj/) 

p0 = R? ' Pb = R?{l + Xyf ■ (20) 

The three cases are studied in the following: 
3.1 Maximum mass for a given radius 

Here maximum compactness implies maximum allowed mass within a given radius of a 
star. For example, if A = 2 and radius 6 = 10 km, the maximum mass is ~ 2.33 M and 
for values of A ~ 100, M max = 2.45 M©. In Table 1, the values of (6 2 / R 2 ) m ax , (M/b) ma x 
and M max for different radii are given for different A. We have plotted the variation of 
maximum compactness {M/b) max and (b 2 /R 2 ) max with A in fig.l and fig.2, respectively. 
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3.2 Maximum mass for a given surface density(p fc ) 

To find the maximum mass, we write the mass in terms of surface density 

(l + A)3/y/2(3 + Ay)V 2 

2^(1 + Ay)2 • W 

To obtain the maximum mass of a realistic star one has to put some physically acceptable 
lower bound on the value of surface density. Rhoades and Ruffini [2] obtained a maximum 
neutron star mass ~ 3.2M & allowing an uncertain nature in the EOS above a fiducial 
density 4.6 x W 14 gm cm' 3 . If we take 4.6 x 10 14 gm cm~ s as the surface density in equation 
(21) we get M max ~ 3.01Mq which is comparable with the result obtained by Rhoades and 
Ruffini [2]. Kalogera and Baym [42] obtained an upper bound on the neutron star mass 
equal to 2.9M Q regarding the EOS as valid upto twice nuclear matter saturation density, 
Pnuc = 2.7 x W 14 gm cm~ 3 . In Table 2 we have shown the maximum possible masses and 
radii for different choices of surface density. 

The EOS obtained by Dey et al [8] for strange stars, namely SSI and SS2, give rather 
low values of maximum mass (M max = 1.33M for SSI and M max = 1.44M Q for SS2 ) 
where the approximated linearised EOS has ^ = a, a constant. For SSI, a = 0.463 and 
pb = 1.15 x 10 15 gm/cm 3 and for SS2, a = 0.455 and pb = 1.33 x 10 15 gm/cm 3 . If we 
consider the same values of pb in equation (21) in our model with (^)o = 1> we get higher 
values of maximum masses and radii as M max = 1.9M Q and b = 7.98 km for SSI and 
M max = 1.77M and b = 7.43 km for SS2, respectively, where A = 53.34 for SSI and 
A = 230.58 for SS2 [28]. A lower value of a will lower the value of maximum mass possible, 
vide Table 4. Note that for a large A, the EOS in our model is almost linear, giving a 
constant value for The variation of M max with pb is shown in fig.3. 

3.3 Maximum mass for a given central density (p ) 

The mass in terms of central density can be written as 

% , Vsa + a) 3 /v /2 

M = y ' — ^— (22) 
2^(1 + Xy) K ' 



The variation of M max with central density is shown in fig.4. Very little information is 

available regarding the central densities of realistic compact stars. However, the central 
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A 


p h = 5.4 x 10 14 


Pb = 10.8 x 10 14 


p h = 4.6 x 10 14 


p h = 5.1 x 10 14 




71 T 


7 

Umax 




7 

®max 


^max 


7 

Umax 


Tl if 

■tVJ-max 


7 

®max 


1 


2.61 


12.19 


1.84 


8.62 


2.83 


13.21 


2.68 


12.55 


2 


2.78 


11.93 


1.96 


8.44 


3.01 


12.93 


2.86 


12.28 


3 


2.78 


11.6b 


1.96 


8.24 


3.01 


12.63 


2.8o 


12.00 


4 


2.76 


11.47 


1.95 


8.11 


2.99 


12.43 


2.84 


11.80 


5 


2.74 


11.33 


1.94 


8.01 


2.97 


12.28 


2.82 


11.66 


7 


2.71 


11.15 


1.92 


7.88 


2.94 


12.08 


2.79 


11.48 


10 


2.68 


11.00 


1.90 


7.78 


2.91 


11.92 


2.76 


11.22 


20 


2.64 


10.79 


1.87 


7.63 


2.86 


11.70 


2.72 


11.11 


50 


2.61 


10.66 


1.85 


7.54 


2.83 


11.55 


2.68 


10.97 


100 


2.60 


10.61 


1.84 


7.50 


2.82 


11.50 


2.67 


10.92 


200 


2.59 


10.59 


1.83 


7.49 


2.81 


11.47 


2.67 


10.90 



Table 2: Maximum mass i^^irnax) in -^^0 corresponding radius bf^ax 

in km of a star for 

different choices of the parameter A and surface density in units of gm cm~^. 



density is useful in determining the static hydrostatic stability of realistic stars under radial 
perturbations. It is known that in the mass-central density diagram, the configurations 
with p > p™ (p™ corresponds to the maximum mass) are unstable with respect to radial 
perturbations. Introducing a time-dependence of the perturbations complicates the cal- 
culations, but the configuration remains unstable with respect to fundamental modes of 
pulsation if the EOS is not changed by perturbation. It has been shown elsewhere ([34], 
[28], etc.) that the Vaidya-Tikekar stars are stable against small radial perturbations. 

4 Discussions 

We have calculated in this paper the maximum mass of a class of compact stars satisfying 
necessary physical constraints. Results obtained are outlined as follows: 

• We have obtained an upper bound on the compactness (M/b) of a star described by 
Vaidya-Tikekar [24] model (see Table 1) for different values of the parameter A. The 
compactness, for large values of A (e.g., A ~ 100), becomes almost independent of A, 
the limiting value being ~ 0.3615. This is less than the Bondi [36] limit ~ 0.3906. 
Negi and Durgapal [37] have shown that in the case of Tolman's type VII solution, a 
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T~> . C , . 

Keierence 


Model 


M max /M Q 


Ruffini and Rhoades [2] 


Neutron star (causality principle) 


3.2 


Haensel et al[15] 


Neutron star (causality principle) 


3.0(5 x 10 14 gm cm' 3 / pb) 1 ' 2 


Baldo et al [40] 


Neutron star(BBB2 EOS(npep)) 


1.92 


Kalogera and Baym [42] 


Neutron star(WFF88 [43] EOS) 


2.9 


Witten [4] 


Strange star( Bag model) 


2.0(B /B) 1 / 2 , B = 56MeV fm~ 3 


Burgio et al [41] 


Strange star(Bag model) 


1.45 - 1.65 


Banerjee et al [19] 


Strange star(Bag model) 


1.54 (5 1 / 4 = 145 MeV) 


Mak and Harko [39] 


Charged strange star (Bag model) 


2.86 


Harko and Cheng [21] 


Strange star(Bag model) 


1.83 


Cheng and Harko [22] 


Strange star (Bag model) 


2.016 


Knutsen[34] 


Vaidya-Tikekar model [24], A = 2 


3.0 


Present work 


Vaidya-Tikekar model [24], any A, 


3.01 (A = 2, p b = 4.6 x 10 14 gm cm' 3 ) 




(J)o = l 


2.82 (A = 100, p b = 4.6 x 10 14 gm cm' 3 ) 




2.45 (A = 100, b = 10 km) 
1.96 (A = 100, b = 8km) 
1.47 (A = 100, b = 6 km) 



Table 3: Maximum mass configurations obtained in different models. 



fluid sphere configuration remains stable for -r- < 0.3428. Based purely on causality 
limit, the upper bound on compactness obtained here, lies slightly above this value. 

• We have obtained a better bound for the radial parameter b 2 /R 2 as compared to the 
earlier results obtained by Mukherjee et al [29]. The new bound obtained here is in 
agreement with the earlier work done by Knutsen [34] for a specific value A = 2. 



• It is observed from Table 2 that for a given surface density (pb), maximum mass first 
increases and then decreases with A. For example, for pb = 4.6 x 10 14 gm cm~ 3 , this 
maxima is obtained when A w 3 and the corresponding values of maximum mass and 
radius are 3.O1M and 12.63 km, respectively. 

• The maximum mass and corresponding radius decreases with the increase of surface 
density A comparative study of maximum masses for different models are listed in 
Table 3. 
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(K) 

\ b /max 


p b = 4.6 x 10 14 


p b = 5.6 x 10 14 


p b = 1.15 x 10 15 


p b = 1.33 x 10 15 






M max 


bmax 


M max 


bmax 




bmax 




bmax 


i 


0.3615 


2.82 


11.50 


2.55 


10.42 


1.78 


7.27 


1.65 


6.76 


0.500 


0.3206 


2.47 


11.38 


2.24 


10.31 


1.56 


7.19 


1.45 


6.69 


0.463 


0.3129 


2.40 


11.34 


2.18 


10.27 


1.52 


7.17 


1.41 


6.67 


0.455 


0.3110 


2.38 


11.33 


2.16 


10.26 


1.51 


7.16 


1.40 


6.66 


0.333 


0.2658 


1.98 


10.99 


1.79 


9.96 


1.25 


6.95 


1.16 


6.46 



Table 4: Maximum mass (M max )m M and corresponding radius b max in km of a star 
for different choices of the parameter a and surface density p b in units of gm cm~ 3 . We 
assumed A = 100. 



If we consider a class of EOSs so that (^)o = o> with, a < 1, the relevant causality 
condition, eqn (17) gets modified to 



A , ^ (1 + A) 



V' 2 °- 2v/A 



v / XTl- ^(20 a + 1)A + 1 



(23) 



It is seen that the maximum compactness, for a given A, decreases as the value 
of a is lowered. For given values of surface density and A, we have calculated the 
maximum mass and the corresponding radius for different choices of a, and these 
are shown in Table 4. It shows that for a star with an EOS giving a = 0.455 and 
having surface density p b = 1.33 x 10 15 gm cm~ s , the maximum possible mass will be 
M max ~ 1.4 Mq and radius b ~ 6.6 km. Although the present model covers only a 
class of compact stars, the results obtained are in conformity with the general results 
obtained directly by numerical calculations. Thus, if one chooses p b = 4 B with the 
Bag constant B = 10 14 gm /cm 3 (56 MeV fm~ 3 ) and a = 0.333 in the Bag model, 
one gets M max « 2M Q and corresponding radius b max rs 11 km (vide [4], [44]). If we 
use the same inputs for p b and a = (^)o an d A = 100, we get M max ps 2.12M q and 
corresponding radius b max ~ 11.78 km, comparable to the numerical values obtained 
by [4], [44]. 



• Consideration of rotation and the presence of charge may change this scenario to some 

extent. However, both the effects will increase the maximum mass limit. For example, 

12 



A 


a 


Umax 


\ b Jmax 


M max (M Q ) 


bmax {km) 


2 


0.0 


0.4234 


0.3439 


2.78 


11.93 




0.4 


0.4474 


0.3529 


2.87 


11.99 


100 


0.0 


0.0252 


0.3615 


2.60 


10.62 




0.4 


0.0252 


0.3616 


2.60 


10.62 


1000 


0.0 


0.0026 


0.3615 


2.59 


10.57 




0.4 


0.0026 


0.3615 


2.59 


10.57 



Table 5: Maximum mass configurations for a charged sphere of surface density pb = 5.4 x 
10 14 gm cm~ 3 . 

Mak and Harko [39] have shown that the maximum mass and radius of a charged 
strange star with a linearized EOS in the MIT bag model are M max = 2.86M and 
b = 9.46 km, respectively, where the value of (dp/dp) = 1/3. In the uncharged limit, 
the corresponding values are M max = 1.9638M@ / y/ Bqq and b = 9.46 km/ 't/Bqq, with 
a scaling with respect to bag constant B = 60 MeV fm^ 3 [4], [15]. This shows 
that in the presence of charge the maximum mass increases which is expected as 
the gravitational attraction is now partly balanced by the Coulomb repulsion. In an 
earlier work [38], the solution obtained by Mukherjee et al [29] was extended to the 
case of a static charged spherical distribution of matter. The solution obtained there 
may be used to study similar effects of charge on the maximum mass. It can be shown 
that for a small A, both maximum mass and the corresponding radius increase in the 
presence of charge, e.g., for pi, = 2p nuc , A = 2 and a = 0.4 (a is a measure of charge 
density [38]), we found that M max = 2.87 M and b = 11.99 km, both of which are 
higher than the values in the uncharged (a = 0) shown in Table 5. However, 

if we consider A = 1000, the maximum mass and corresponding radius for the same 
values of pb and a are M max = 2.59 M Q and b = 10.57 km, exactly the same as in 
the uncharged case. Thus, the effect of charge is very small on maximum mass and it 
becomes negligible for higher values of A in this model. Similarly, if rotation is taken 
into account, maximum mass limit is expected to increase marginally as pointed out 
in ref[15]. 
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To conclude, we have shown that the Vaidya-Tikekar model and its solution of Mukherjee 
et al [29] provide a simple method of studying systematically the maximum mass problems 
of compact stars. The calculations can be done without any prior knowledge of the EOS. 
Although the class of stars described by this simple model may not cover, in principle, 
all compact stars, it nevertheless captures the gross features of known compact stars, as 
is evident from a comparison with results obtained by solving the TOV equation with 
different conventional equations of state for Neutron or strange stars, given in Table 3 and 
Table 4. 
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Figure 1: Variation of (M/b) max with A (solid line for a = 1 and dashed line for a = 0.455). 



Figure 2: Variation of (b 2 /R 2 ) max with A (solid line for a = 1 and dashed line for a = 0.455). 
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Figure 3: Variation M max /M & with surface density pb (gm/cm 3 ) (solid line for A = 3 and 
dashed line for A = 100). 



Figure 4: Variation M max /M & with central density po(gm/cm 3 ) (solid line for A = 3 and 
dashed line for A = 100). 
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